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Abstract. A general treatment of vorticity-free, perfect fluid perturbations of Kantowski- 
Sachs models with a positive cosmological constant are considered within the framework 
of the l-|-l-|-2 covariant decomposition of spacetime. The dynamics is encompassed in six 
evolution equations for six harmonic coefficients, describing gravito-magnetic, kinematic and 
matter perturbations, while a set of algebraic expressions determine the rest of the variables. 
The six equations further decouple into a set of four equations sourced by the perfect fluid, 
representing forced oscillations and two uncoupled damped oscillator equations. The two 
gravitational degrees of freedom are represented by pairs of gravito-magnetic perturbations. 
In contrast with the Friedmann case one of them is coupled to the matter density pertur¬ 
bations, becoming decoupled only in the geometrical optics limit. In this approximation, 
the even and odd tensorial perturbations of the Weyl tensor evolve as gravitational waves 
on the anisotropic Kantowski-Sachs background, while the modes describing the shear and 
the matter density gradient are out of phase dephased by 7r/2 and share the same speed of 
sound. 
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1 Introduction 


The observed large scale distribution of galaxies, the fluctuations about the isotropic cosmic 
microwave background radiation and the late time acceleration of the universe seems to be 
well described by the ACDM model, which is based on the assumption that the geometry 
of the universe is given by the Robertson-Walker metric - see e.g., [1-5]. However this fit is 
not perfect [6-9] and because more than 95% of the matter budget needs to be described by 
the dark sector, it is worth exploring what effect alternative cosmological models have on the 
basic properties of the Universe [10-21], 

It is known that anisotropies in the Hubble and deceleration parameters cannot be ex¬ 
cluded by present observations [22], [23] and [24]. Consequently, perturbations of anisotropic 
cosmological models have been considered by many authors, e.g., [25-32], using both gauge 
dependent methods (e.g., [33]) or Bardeen’s gauge invariant formalism [34]. For example the 
perturbations of homogeneous and anisotropic universe of the Bianchi I type was investigated 
in [31, 32] by using Bardeen’s gauge-invariant method. However, the variables in Bardeen’s 
theory are defined with respect to a particular coordinate system, making their geometrical 
and physical meaning not very transparent See the discussion in [35]. By using a covariant 
approach, one circumvents these problems by using the spatial curvature rather than the 
metric as the defining variables [36, 37]. In this way, a set of gauge-invariant perturbation 
variables can be easily identified as the ones that vanish on the chosen background [38- 
44, 50]. This feature of the covariant approach makes it a very versatile method for studying 
perturbations on a variety of backgrounds and physical situations and relating the results 
obtained in a unified way [45-49] 

In this paper we present for the first time a general treatment of the vorticity-free 
perturbations of Kantowski-Sachs cosmologies with positive cosmological constant, extending 
earlier work [51], which focused only on the scalar perturbation sector. Here we present for 
the first time an analysis of a full scalar, vectorial and tensorial perturbations, focusing on 
gravitational and matter wave evolutions. 

In order to achieve this, we use a covariant and gauge invariant method [52], in which 
spacetime is first split into a 1-1-3 form. The formalism has been mainly employed for com¬ 
putation of cosmological perturbations on a Friedmann background, applying the standard 
decomposition theorems [35] (see for example Refs. [38, 40, 53, 54]). If the 3-space at 
each point has a unique preferred direction, a further decomposition of the spacetime into 
a l-|-(l-|-2) form is useful in situations where spacetime admits a spherical or Locally Ro¬ 
tational (LRS) symmetry. This was first employed in Ref. [55], where the spatial direction 
was singled out by local rotational symmetry (LRS). The formalism was developed with the 
purpose of investigating general gauge-invariant perturbations of the vacuum Schwarzschild 
spacetime [56]. With the further generalisation presented in Refs. [58] and [59], it became 
possible to describe gauge-invariant perturbations of LRS class H spacetimes for which the 
complete set of evolution and constraint equations are given in Ref. [60]. 

In this paper, the variables describing an almost Kantowski-Sachs spacetime are ex¬ 
panded into harmonics. We find that the perturbation dynamics is described by six evolution 
equations for six harmonic coefficients, together with a set of algebraic expressions, which 
determine the evolution of the rest of the variables. The evolution equations can be split into 
two sets - four which are sourced by the perfect fluid, representing forced oscillations and 
the remaining two describing damped oscillating gravito-magnetic perturbations. We fur¬ 
ther analyse the equations using the geometrical optics approximation and find that four of 
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the gravito-magnetic quantities evolve as gravitational waves propagating on the anisotropic 
Kantowski-Sachs background, while the shear and the matter density gradients, which are 
out of phase by 7r/2, share the same speed of sound. 

The paper is organised as follows. In Section 2 we briefly review the 1+3 and 1+1+2 
covariant approaches. In Section 3 a Kantowski-Sachs type background hlled with a perfect 
fluid is introduced. The equations governing the linear perturbations are derived in Section 
4. All type of perturbations (scalar, vector and tensor) are investigated, however they are 
restricted by vanishing anisotropic pressure and energy flux, which mean that we assume a 
perfect fluid and that the 1+3 split is done with respect to the 4-velocity also in the perturbed 
spacetime. For simplicity we also choose to put the vorticity to zero. This imply that the 
hypersurfaces perpendicular to the 4-velocity are well defined. The gauge degrees of freedom 
in the choice of frame are analysed in Appendix B. Applying the commutation relations given 
in Appendix C and the useful relations for the vector and tensor spherical harmonics given 
in Appendix D, the equations governing the perturbed system are derived in Appendix E. 
Then by fixing a frame we find that the perturbed spacetime can be described by six type 
of harmonic coefficients. The evolution equations for these six variables are given in Section 
4. The behaviour of the perturbations in a geometrical optics approximation is discussed in 
Section 5, while Section 6 contains some concluding remarks. 

Units where SvrG = 1 and c = 1 are used throughout this paper. 

2 The 1+3 and 1+1+2 covariant formalisms 
2.1 The 1+3 covariant formalism 

Let be a time-like vector field obeying the usual normalisation condition u°‘Ua = — 1 and 
hab a spatial 3-metric satisfying u^'hab = 0. Then the 4-metric gab can be decomposed as 

9ab UaUb T hab • (^d) 

We denote the 4-dimensional (4D) and 3-dimensional (3D) volume elements by rjabcd = 
^O'bc = Vdabcu'^, respectively. Angular brackets ( ) on indices denote 
symmetrised and trace-free tensors which are projected in all indices with the metric hab- 
Round brackets ( ) and square brackets [ ] on indices denote the symmetric and antisymmet¬ 
ric parts, respectively. A dot denotes covariant derivatives along the integral curves of tt“, 
while Da is the projected spatial derivative 

Tb..c = u°'VaTb„c , ( 2 . 2 ) 

DaTb..c = ha%\.ha^Vdn.j . (2.3) 

For vanishing vorticity of u“. Da is the 3D covariant derivative compatible with the metric 

hab- 

The kinematic quantities are introduced through the decomposition of the 4D covariant 
derivative of as ^ 

^a'^b ^ab T ~^^hab + ^ab '^a-^b ; (^d) 

where aab = D(^aUb) is the shear, 0 = D“ua the expansion and Uab = the vorticity of 

tt“, finally +„ = U(^a) = ha ^Ub is its acceleration. Since Uab is space-like and antisymmetric, 
containing 3 independent components, we introduce its (Hodge-) dual u!a = ea'^bJbcl‘^- 
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The gravito-electro-magnetic quantities arise from the 1+3 covariant decomposition of 
the 4D Weyl tensor Cabcd as 


Eab = CacbdE'u^ and Hab = ""^CcdbeU^- (2.5) 

The quantities Eab and Hab are the magnetic and electric parts of Cabcd, respectively. The 
Weyl tensor is given^ [36] by 

'^Cabcd ^c'^[aEb]d '^d'^\aEb\c T .^c[a^b]d Ed,\ahb]c ^ab Hj^^aUb] ■ (2.6) 

As well-known, the Weyl tensor (posessing the symmetries of the Riemann tensor with 20 
independent components, and also being traceless, meaning 10 conditions) has 10 independent 
components. Its electric and magnetic parts have each 5 independent components, as they 
are expressed by traceless and symmetric 3-tensors. 

The energy-momentum tensor Tab of matter fields is decomposed with respect to an 
observer with 4-velocity in the standard way: 


Tab = fJ-UaUb + 2q^aUb) + PKb + T^ab ■ 


(2.7) 


The quantities fi, qa, p and TTab are the energy density, the energy current vector, the isotropic 
pressure and the symmetric, trace-free anisotropic pressure tensor of matter. 

The Riemann tensor can be expressed by the metric components, gravito-electro-magnetic 
quantities and matter variables as follow. First we use the decomposition of Rabcd into its 
Weyl and Ricci {Rab) contributions 

R 

Rabcd Cabcd T 9a[cRd[b 9b[cTd]a ~^9a[c9d]b , (2.8) 

where R is the Ricci scalar. Then by applying the Einstein equation 

T 

Rab ~ Egab + Tab ^5af) ; (2.9) 

with T = g'^^Tab and cosmological constant A, and Eqs. (2.6)-(2.7), we find 


Rabcd 2 (^UcUi^aEbjd '^d1^[aEb]c T Ec[a^b]d Ed[a^b]c ^ab T[i[c^d] ^cd Ri[a^b]} 


(2A — p — 3p) 


(uaU^c^d]b ^6^[c^d]a) T 


2 (A + /i) 


ha[c^d]b 


^ \ tij L/ L< CiJ Ci _/ ' ^ 

+9a^[c^d]fe QbU^chd]a T ^bQ[chd]a T 9a[c^d]b 9b[c^d]a • 


( 2 . 10 ) 


The full set of equations arise from the Ricci identities for and from the 4D Bianchi 
identities and can be found for instance in Refs. [52], [61]. 

A 3D curvature tensor can be defined in the following way (see [62]): 

1 (3) . 

2 RabcdV^ = D[aDb]Vc - ^^abV(c) , ( 2 . 11 ) 

where is an arbitrary 3-vector. In the vorticity-free case (as a consequence of Frobenius’s 
theorem) Rabcd is the Riemann curvature of the hypersurface with metric hab- Alternatively 
it can be given in terms of the Gauss’ equation [63]: 

Rabcd = hahjh^hd Rabcd “ {DcUa){DdUb) + {DdUa){DcUb) - (2.12) 

^The definition of Eab differs by a sign in [36]. 


_ 4 „ 





From Eqs. (2.12), (2.4) and (2.9) the curvature scalar Rabcd is [62]: 

i^)R 02 

-—= ^, + A- — + a-^aab - ooaboo^^^ (2.13) 

giving the usual Ricci scalar on a spatial hypersurface orthogonal to (when ujab = 0). 

2.2 The 1+1+2 covariant formalism 

When there is a unique preferred spatial direction at each point, it is worthwhile performing 
a further decomposition of spacetime into its so called 1+1+2 form. The preferred spatial 
direction will be singled out by a normalised vector held n“ (n“na = 1 , n°‘Ua = 0) and this 
allows us to decompose the metric hab as 


hab — fT'a'fT'b + ^ab j (^'1^) 

where N^b is the induced 2-metric on the surface perpendicular to both and The 
alternating Levi-Civita 2-tensor is dehned as £ab = ^abc'n-’^- Curly brackets { } on indices will 
denote symmetrised and trace-free tensors which are projected in all indices with the metric 
Nab- A bar on vector indices will denote projection onto the 2-sphere: Va = NabV^. The 3D 
covariant derivative can then be projected into two parts as [60] 

fb..c = n^DaTb..c , (2.15) 

5aTb..c = Na^^N^K.Na WdR,,, . (2.16) 

A further decomposition of the 1+3 vector and tensor variables with respect to n“ are given 
according to Ref. [60]. The kinematical vectors and tensor can be decomposed as 

+“ = ^n“ + .A“, (2.17) 

a;“ = Dn“ + D“, (2.18) 

aab = E (^UaUb - ^^ab^ + 2S(„n6) + Sab , (2.19) 

while the gravito-electro-magnetic variables are 

Eab = £ (^riaUb - ^^ab^ + ‘2-£(anb) + £ab , ( 2 . 20 ) 

Hab = H ^UaUb - + ‘NhLf^a'n-b) + Hab ■ ( 2 - 21 ) 

In the above decompositions the scalars, the 2-vectors and the symmetric 2-tensors represent 
one, two and two independent components each, respectively. Finally, the 1+2 form of the 
energy current vector and anisotropic pressure tensor are 


Qa — Qna + Qa i ( 2 . 22 ) 

T^ab = n (^riaUb - ^Hab^ + 2n(anb) + Flab ■ (2.23) 

Here S“, "H®, <5“ and n“ are 2-vectors and S^b, £ab: Hab and H^b are trace-free, 

symmetric 2-tensors perpendicular to both and n“. 
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The additional fundamental variables of the 1+1+2 formalism arise from the projected 
time derivative of n“ and the 3D covariant derivative of Ha- The time derivative can be 
written as ha = Bua + aa, where Ua is a 2-vector (since han“ = 0). Now, from u°'na = 0, it 
follows that = —naA°'. Using (2.17) one obtains B = A. Hence 

ha — AUa + Ola (2.24) 

{oa = ha)- The 3D covariant derivative of Ua can be decomposed as: 


DaTlb 


Cab T 2 ^ab T C^ab T 


(2.25) 


where (f) is the sheet expansion, = n’^DcUa is the acceleration, ^ab is the shear of n“ and 
C represents its rotation in the local 3D space. Here Ua and Oq are 2-vectors and Cab is a 
trace-free, symmetric 2-tensor perpendicular to both and n“. 

The full set of evolution and constraint equations for perturbed LRS spacetimes are given 
in Ref. [60]. There are no evolution equations for A, Aa, oia and there is no propagation 
equation for Oq. These are determined by fixing a particular frame [60]. 

We define the 2-dimensional (2D) curvature tensor TZabcd as 


^TlabcdV^ 


^[a^b]^c ^^ab^c C^ab^c , 


(2.26) 


where U“ is an arbitrary 2-vector. Similar definitions are given in [64] for higher dimensional 
spacetimes. For vanishing H and C it agrees with the usual Riemann curvature tensor of Nab- 
This definition gives (see Appendix A) 

TZabcd = NlNlN^N\Rijki+{5aUd) {5hUc)-{5aUc) {6bUd)-{Sand) {Sbnc)+{5anc) {hud) , (2.27) 


where Riju is the usual 4D Riemann tensor. By using Eqs. (2.10), (2.4) and (2.17)-(2.25), 
TZabcd is expressed as 


TZabcd — 


1 


-n- 2 ^--s-— + 


2 (A + /i) 


+-^cfa S + 2^6+ + 1 ^ “ 


20 V 


3 / 


^claA^b]d 


b]d 


20 


20 


^b]d + 4>Cb]d + 


~^d[a + ‘2£b]c + ^ j ^b]c + ‘t’C^c + 

+2ec[a [£^^b]d - '^Cb]d] - 2erf[a [HS^jc - CCb]c] 

+ 2 (n^ — V) ^c[a^b]d + 2 (Cc[aCb]d “ • 

Defining the 2D curvature scalar TZ = TZabcd, we hnd 

20 \ 




^b]d 

^b]c 


(2.28) 


7^=-(^ + A)-^-2^--(S-—J +^-2(H2-V)+S,6S“'-Ca6C' • 
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(2.29) 









3 The Kantowski-Sachs background 

The spatial sections of Kantowski-Sachs cosmologies have topology R x S'^ and are the only 
spatially homogeneous cosmologies that do not fit into the Bianchi classification. This is 
due to that their isometry group does not admit a 3-dimensional subgroup that acts simply 
transitive on the hypersurfaces of homogenity. These metrics are Locally Rotationally Sym¬ 
metric (LRS) and belong to the LRS class II, characterised by coab = ^ = Hah = 0 [55, 65]. 
The square of the line-element can be written as 

ds^ =—dt^ + a\{t) dz^ + a^{t) + siv?-ddip^) . (3.1) 

The 4-velocity of comoving observers is u = d/dt and the direction of anisotropy is n = 
a^^d/dz, where z is dimensionless. The coordinates d and tp are the polar and azimuthal 
angles on respectively. The scales oi and 02 have the dimension of time, 02 being assumed 
to be sufficiently large to avoid periodic structures in the angles emerging in the observable 
Universe. Symmetry and normalisation implies [55]: 

hfe = n^-DaUb = 0 , hj = 0 , (3.2) 

i.e., Ua is geodesic on local 3-space with metric hab and is Fermi propagated along the integral 
curves of In the spacetime given by Eq. (3.1) the non-vanishing kinematical variables of 
formalism are the expansion [51]: 

e = L + 2^ , 

Oi 02 

and the scalar part 

3 V«i «2/ 

of the shear aab- Given an equation of state p = p (p) for the pressure p and energy density 
p of the perfect fluid, for any given cosmological constant, the Kantowski-Sachs models are 
completely determined in terms of the shear S, expansion 0 and p. The electric part of 4D 


Weyl tensor is then determined algebraically as (see Eq. (100) in [60]) 

3T = -2(/r + A)-3S2 + ^02 + S0 . (3.5) 

O 

The evolutions of S, 0 and p are governed by Eqs. (96), (94) and (95) of [60]: 

p = -Q{p + p) , (3.6) 

0 = -^-^S2-l(^ + 3p)+A, (3.7) 

t, = - {p + A) + — S0 — -0^ , (3.8) 

where we have used Eq. (3.5). 

For Kantowski-Sachs spacetime the 2D scalar curvature (2.29) becomes 

K = ?(, + A)-2f-i(E-f)T2(, + A) + fE^-?|f = ^. (3.9) 


(3.3) 

(3.4) 
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that is two times the Gaussian curvature of the 2-spheres. Taking the time derivative of IZ, 
and using Eqs. (3.6)-(3.8), we find 

(3.10) 

One of the evolution equations (3.6)-(3.8) can be replaced by (3.10). 

In summary the non-vanishing quantities on the background are given by the set 

= (3.11) 

or equivalently 

5® = {0,S,7^,/r,p}. (3.12) 

To zeroth order (on the background) £ (or TZ) are given in terms of the other quantities. 

General orthogonal spatially homogeneous LRS class II spacetimes emerge as slight 
modifications [65]. The Kantowski-Sachs cosmologies are the only ones with TZ > 0. If 7?. < 0 
the spacetimes are of Bianchi type III and the only modifications to the above equations are 
to replace sini? by sinht? in equation (3.1) and to change TZ to TZ = —2la\. For TZ = 0 there 
are solutions of Bianchi type I/VIIq. Due to that (3.9) now becomes a (satisfied) constraint, 
one of the evolution equations can be dropped. There are also TZ = 0 with S = T = 0. For 
these the sheet expansion cj) is in general nonzero and the system is given by equations (3.6) 
and (3.7) plus the constraint 

//+ A —-0^= 0 . (3.13) 

These are flat (when c;/!) = 0) or negatively curved Friedmann models, which also fall into the 
Bianchi I and V classes respectively. 

In this paper we consider perturbations of the TZ > 0 models, but the above discussion 
shows that perturbations of the TZ < 0 can be done in an analogous way. Perturbations of 
Friedmann models and Bianchi I models have been considered elsewhere [38-44]. 

4 Vorticity-free, perfect fluid perturbations of Kantowski-Sachs cosmolo¬ 
gies 

For simplicity, we assume the perturbed fluid is irrotational and we use a frame associated 
with the fluid. This requirement assigns the reference 4-velocity in the perturbed spacetime 
and in a 1-1-3 covariant formalism the frame is completely fixed. In this frame the energy 
current qa and ui^ = vanish. Moreover, we neglect the anisotropic pressure contri¬ 

butions to the energy-momentum tensor, i.e., iTab = 0, restricting our analysis to barotropic 
perfect fluids. 

The variables of 1-I-1-I-2 formalism are defined with respect to the frame vectors u and 
n. Therefore the variables are not frame-invariant in general (see Appendix B). (Of course 
their combinations may result in frame-invariant quantities [56].) The frame choice does not 
fix completely the mapping between the perturbed and the background geometry [38], [66], 
[67]. The variables vanishing on the background are invariant for the remaining gauge fixing 
in this map according to the Stewart-Walker lemma [50]. Therefore instead of T, 0, S, fj, 
and p we use their gradients 

Xa = 6a£ , Pa = 5aD , Wa = SaO , 

fJ-a = SaP , Pa = ^aP , (4.1) 
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that vanish on the background. As will be shown in section 4.2, the hat-derivatives (2.15) 
are related to the da derivatives when the vorticity vanishes. 

Hence we have the following nonzero first order quantities (that vanish on the back¬ 
ground) : 

= {Xa,Va,Wa )/^a; Pa;-A, vAa, C) Cafe} • (^■^) 


From the l-|-l-|-2 equations in [60] the following nontrivial evolution equations then hold 
on the perturbed Kantowski-Sachs spacetime 




— .A ) -|- S°'aa , 


3/12 


2i= + + n , 


20 \ 




3 J 


1 


20 \ 


fla = -{T, - —] + Pa) - ilJ'+P)Wa + flAa , 


3 / 




3£ 


(4.3) 

(4.4) 

(4.5) 

(4.6) 


X-a = 2l S- ^ ) Xa + Y[^-- f (/^a + Pa) + £Aa + £^‘^Ubnc , (4.7) 


Va-^Wa= ^ ( S- 


20 


1 


20^ 


Va--Wa -Xa + - (pa + 3pa) + S- — A " U^Ab, (4.8) 


^{ab} ^ ^ ^ab 4“ £ab j 

C{ab} 2 2 ^ Cab 4 “ *^{a®b} £c{aT^b} ' 

The equations containing both propagation and evolution contributions are 

Wa-5aA= (^-Q]Wa- 3SHa - ^ (pa + 3pa) + ©A + Sad’^Ab 


Q^a da — ( ^ 4“ 2 ^ {Aa 4“ fla) £abT~L j 
40 \ 

2Sa — Aa = da A — ( S -|- —j — 3Saa — 2Ta , 


1 


^a+TT^ab^ =T S- — Ta + 


40' 


{3£ — 2p—2p) 


3£ 


1 


40 \ 


3£ 


3£ 


1 


c-b cc 


n-a- -eab£^ = - ( S - — j 7/a - ^^ab^'’ - —£abA^+ -^£aba^ - ^ebcd^£ 


(4.9) 

(4.10) 

(4.11) 

(4.12) 

(4.13) 


Sa + (4-14) 


(4.15) 


^There are some minor misprints in [60]. In Eq. (36) the term (Ea — tab)^^dj) should probably read 
—2ea6ll*’b, in Eq. (40) (Ea - eac)II°'i/’b —t —2eaAl‘'7’b, in Eq. (52) the term —(§0 -t |E)Ea6 —>■ —(§0 - E)Ea6 


-t (§0 -t E)(^a — aa) —>■ ■p^tabOi’ + (^ 

-\T.eabH*’ and in Eq. (80) -(§0 - iE)(E, - tab^*’) ^ +^0 - - eab^i’’). 


in Eq. (53) the terms —^eabO. + (^0-|-E)(^a — fla) —S- -l-^£a6a + (|0-t E)(^a -t fla), in Eq. (76) —'Seab'H 
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^{ab} ^c{a^b} 2 ^ 


20 _ rc-^f (3f + /i+p)^ 

^ab £c{a^ ’^b} n ^abi 


^{a6} “1“ ^c{a^b} 


20 


3£ 


H ^ ) T~l-ab + -;r£c{aCb} + ^c{a^^Sb} ■ 


The pure propagation equations are 


20 


0 


= - S-— S + - + 5“aa -£- 


2 (/r + A) 


2f = , 

n = -5^Ha , 

*^a — 5 

Pa = - ifJ'-\-P) SaA , 
y - 

9,^_, ^ 35] 


1 


3S 


Sa = 77 K + TT^a - , 


2£a = Xa- 3£aa - 3£:eab-H^ - 25^£ab+ -/^a , 

in-a = 5aU - 3£eab£‘^ + 3Se,fef^ - 26^Hab , 

^ 3S 

^{a6} *^{a^6} “I TT^ab ^c{aT~^b} ’ 


C{a 6 } ( ^ 2 ) *^{ 0 ^ 6 } ^afc 


Finally, the constraints are 


e“^5„A = 0 , 

Pa = - + P) Aa , 

s'^^da^b = -3Se + H , 


ea^n ^Tafe-- 


aE„+f„ 


3 ) 


Va --Wa]+ = -2eabn^ 


(4.16) 

(4.17) 


(4.18) 

(4.19) 

(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 

(4.27) 

(4.28) 

(4.29) 

(4.30) 

(4.31) 

(4.32) 

(4.33) 

(4.34) 


The equations (4.3)-(4.34) were derived from the generic 1+1+2 equations given in Ref. 
[60], by use of the commutation relations given in Appendix C. 
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4.1 Harmonic expansion 

Following Ref. [51] we expand the scalar perturbation variables into harmonics as 


(4.35) 




The coefficients depend solely of time. The function is the eigenfunction of the 

Laplacian A = n“Van^Vb and it is constant on the 2 ; = const hypersurfaces: 

^2 

AP^ll =-, (5aP^ll = P''" = 0 . (4.36) 

af 

Here are the constant comoving wave numbers in the direction of anisotropy and the scale 
factor ai in this direction obeys 

- = S + I . (4.37) 

ai 3 

The harmonics are introduced on the 2-sphere as 




l{l 1 ) 


Q 


.,m 


Ql,m ^ Ql,m ^ Q 


where 6'^ 


5a5°‘, and the second scale factor 02 satisfies 


02 _ 1 
02 2 



(4.38) 


(4.39) 


For a given I value the index m runs from —I to I. Due to the symmetries of the background 
spacetime the index m never appear explicitly, therefore we will use the following notation: 

jagfcx = — ±Qki_ ^ gfcx = = Q , (4.40) 

«2 

with k‘j_ = l{l + 1) comoving wave numbers in the perpendicular direction to n“. 

The vectors and tensors can be also expanded in harmonics by introducing the vector 
and tensor spherical harmonics [56, 57, 59]. The even (electric) and odd (magnetic) parity 
vector harmonics are 

= a 26 aQ'^^ , = a2eab6^Q'^^ , (4.41) 

and the vector 4'^ can be expanded as 

'ka= E ■ (4.42) 

fc||,fcx 

Similarly, the even and odd tensor spherical harmonics are 

Qlt = , Qlt = , (4.43) 

and the tensor 'ifat can be expanded as 

E^'" +%^k^Q':t) ■ (4.44) 

k\\,k± 
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This decomposition of the vectors and tensors encompasses both an expansion into spherical 
harmonics and into even and odd modes, similarly to the Regge-Wheeler decomposition of the 
perturbations of spherically symmetric spacetimes, leading to the Regge-Wheeler equation for 
the odd modes [68] and the Zerilli equation for the even modes [69, 70]. In our decomposition 
however the coefficients exhibit a {t,z) dependence, rather than (t, r). Some useful relations 
involving the vector and tensor spherical harmonics are enlisted in Appendix D. 


4.2 Relations between harmonic coefficients from commutation rules of covari¬ 
ant derivatives 

A Kantowski-Sachs spacetime filled with perfect fluid is characterised by the time-dependent 
scalars G = {/r, T, S, 0,p}. We expand their d^-derivatives into harmonics on the perturbed 
spacetime. Thus, applying the commutation relation Eq. (C.4) for zero-order scalars and 
using (D.6), we find in the absence of vorticities the following relations: 

Pfc||fcx ^ ■ (4.45) 

However, we could also expand the anisotropic direction derivatives of G into scalar 
harmonics, as they are also first-order. This expansion is 


fcyfcx 


(4.46) 


Then the even parity part of Eq. (C.3) gives 




02 


Ol 


fciifc I 


(4.47) 


a constraint on and emerging in the absence of vorticities. Using Eq. (4.45), 

the odd parity part of Eq. (C.3) becomes trivial. 


4.3 Pull set of evolution and constraint equations for the harmonic coefficients 

From the commutation rules we have found that some harmonic coefficients vanish in the 
absence of vorticities (see Eq. (4.45)). There is a further coefficient = 0, the vanishing 

of which follows from Eq. (4.30). Using these relations, the perturbation equations (4.3)- 
(4.34) can be expanded into harmonics and are given in Appendix E. Some of these equations 
are first integrals of the rest. We have found 17 independent constraints for 28 variables, 
as presented in Appendix E. We consider adiabatic matter perturbations p = p (//) giving 
~ the square of the matter speed of sound.^ 

In the frame associated to the fluid {qa = 0) the vorticity oJa also vanishes because the 
fluid is irrotational. Moreover we have assumed -Kab is negligible. These assumptions fixes 
completely the frame in a 1-1-3 covariant formalism, however do not in a 1-I-1-I-2 description 
where a dyad (u“, n“) must be assigned (see Appendix B). The quantities Qa, oJa and -Kah 
are invariant for some part of the infinitesimal transformations which fixes n“. In particular, 
they are invariant under the infinitesimal translations given by la- We have 2 gauge degrees 

^The adiabatic assumption for the total fluid perturbation, leading to = p//i is a good approximation 
when one of the matter components dominates: for instance, = 1/3 in the radiation dominated area and 
Cs ss 0 for the dust dominated regime. 
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of freedom to sign n“ perpendicularly to on the perturbed spacetime and to fix completely 
the frame. In the frame = 0 Eqs. (E.l) and (E.24) become constraints, indicating that 
the fixing of la reduces the degrees of freedom by four. Thus, we have 6 degrees of freedom 
describing fully the vorticity-free perturbations in the adiabatic case with Qa = 0 = T^ab- 
These variables can be chosen as They are 

invariant for the /^-infinitesimal translation. Their evolutions are governed by two sets of 
decoupled equations, which follow from the evolution equations (E.20), (E.28), (E.30), (E.7), 
(E.31) and (E.8) of Appendix E by employing the constraints. 


4.3.1 Uncoupled evolutions of gravitational perturbations 


The two coefficients and form a decoupled system 

—T 3 _rp iku 

Tfc||A:j^= — 2 + — (1 — D) 




with the coefficients 


'2^ 

2ik\\^ I a\ 


ai 


ai 


-^-bc + 9j:e\£ 1. - {2E + F)n^ 


k\\k\ > 


2kl 


fcl 9S 


2kl 


B = —;5- 1 -^ -|--- h 3£ — —;5- 


2 - fc' 




+ 3S ( S + ? 


Qj^ Qjr) 


CLn 


C = B-^ 


2-k] 


+ 3£ 


D = C + 


k- + P 
B 


E=- C- 


£_ 

B 


+ 


W 




(4.48) 

(4.49) 

(4.50) 

(4.51) 

(4.52) 

(4.53) 

(4.54) 


The second equality in (4.50) follows from Eqs. (3.5) and (3.9). 

Equivalently, the system can be rewritten as decoupled second-order linear homogeneous 
ordinary differential equations: 


where 


—T _j. 

'^kuk, +QniBk..kt +QHo'Hk..ki = 0 , 


^.=|(2F + 2F+E«)-4,„i^ 


(4.55) 

(4.56) 

(4.57) 




1-D 

ai 


2kl 


ai 


+ ai (9SE - BC) 


+ 3-(f+ED) 


-3{E+ED) 


|lnL_^-^(2E + F) 

at ai 2 


(4.58) 
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(4.59) 


Q 1 

q^^ = -{2E + 2F + T.Dy-\n 


'2/c^ 

^ + ai iyEE-BC) 


1- D 

2A:n 

d 

l2kf. 1 

2gH0 — 

^ + Ol (9FE - BC) 

- 3 (2E + F) — In 

^ + ai (9SF - BC) 

Ol 

Ol 

dt 

Ol 


+ '^{E+ED){2E + F) + ?,j^{2E + F) . (4.60) 

The equations (4.55)-(4.56) represent wave equations with friction. As will be shown 
in Section 5 of the paper, the quantities obeying these equations represent the gravitational 
wave degrees of freedom. 


4.3.2 Evolutions with matter sources 

The coefficients and also form a system of differential equations coupled 

to the density gradient as follows 


i-y = 

f^k\\k I 


1-3 


B 


~3~ 


r^k\\k I 


^2 


{p + p) 


X 






(4.61) 




k\\k I 


■“kwh I 


02 [P+P) 


Pkuk I ' \ ^ 2 ) ^knk, 




k\\k \ ’ 


3S 


2a2B 


V P+P yT _^(p I pTl'^ 


Pku k 


■ T 
^fcii k 


ik\\ 


-11^ 

) Pkn k 


—T 

n 


Here we have introduced the additional notations 


k\iki iki\ r\pT 


2kl 

" (1_C)- 


2 — 

a| a^B 


5-1 = 


FB 


£ + 


3S2\ 2 - A:^ 


2a| 


-£- 


kl 


3 ^ £ 

+ 2^”^ S 


(4.62) 

(4.63) 

(4.64) 

(4.65) 

(4.66) 


This system is equivalent to the ones describing scalar perturbations studied in [51]. 
Here the variables 


Va = a^, Za = aDaQ, Ta = aDaa\ (4.67) 

fi 

Sa = Da (4.68) 

where a is the average scale factor defined through 0 = 3d/a and Sab is the traceless part of 
the 3-Ricci tensor, were used. When projected onto the 2-sphere and expressed in terms of 
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the variables Ha, 14, Wa and Xa (see (4.1)) they read 


Va = a^, Za = aWa , Ta = \aT.Va (4.69) 

5a = y SXa + a 05 - 0S + K - . (4.70) 


Due to equation (4.45) we only have to consider the even parity components, 

and of the variables (4.1). The three latter are solved for in terms of 

£k k^i ™ equations (E.51), (E.52) and (E.55) in appendix E. Substitution 

of these into equations (C.1)-(C.4) in appendix C of [51] reproduces the system (4.61)-(4.64). 

We proceed with transforming Eqs. (4.62)-(4.64) into second order oscillator equations 
for each of the gravitational perturbations and T-Lk^^k^t with source terms given 

by the matter perturbations induce forced oscillations. In turn, then 

these gravitational perturbations act as sources for the first order evolutions of Eor 

the gravitational perturbations we obtain 




± ' T T '^SO \/ 


(y ' (/ (Z ' ^ ^ A. y/ y 


and 


• T 


ihSy^Q y 


'^kffkj^ I'Hl'^kffk^ 4" 0x027? 


with the coefficients 


(?si=2S + ^ 


(1 - C) B - (fi + p) d 


20 


0 


20 


Sso = ^ - I - (0 + 3S) - ^ In [o 2 {p + p)] 


(4.71) 

(4.72) 

(4.73) 

(4.74) 

(4.75) 

(4.76) 


qsi = Wi l+0 + ^S(l + C)-^ln(oiP) 
—W 2 + — (1 + 3(7) — — In (/i + p) 


(4.77) 






-VE 2 


3S (1 - C) 


d 


+ 1^2 (0 + 3E) 
1“ 


^ 20 d , , , 

^-T + S‘“ 




(4.78) 
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= ^ _ y _ 1.11 Qi-P 

2B {1 — C)B — {fj, + p) \ 3 S dt p + p 


_ d / 3S \ P 
dt\2a2B) 2a2B 


+ 


+ 


W 3 -W 4 


2o2 02 1(1 - C)B-(p+p)j ’ 


and 


~ s ^ ' 

Q , f. ^ f^ + p\ P d 1 

9wo = 5 + ( 1 - 1 - + -- 


0 7S d 


B 


dtS 


y ^ [«2^ (1 - C)] > 


where we have denoted 


VFi = 

W 2 = 


(l-C)B 


(1 - C)B- (p+p) 
(P + P) 

(1 - C)B- (p+p) 


W3 = 




40 S (4 - 3C) 3S (p + p) 


2B 


W 4 = 


4ln(^ + p)+(E-^) 


40 S 

T~J 




In the next section we analyse the high frequency limit of these equations. 


(4.79) 

(4.80) 

(4.81) 

(4.82) 

(4.83) 

(4.84) 

(4.85) 

(4.86) 

(4.87) 

(4.88) 


5 Geometrical optics approximation 


In this section we follow Isaacson’s definition [71], [72] of gravitational waves on a curved 
background in a geometrical optics approximation. The key concept is that gravitational 
waves are periodic perturbations with a wavelength much shorter than the curvature radius 
of the background. This is known as the geometrical optics approximation, or the high 
frequency limit. In the notations of the present paper the physical wave numbers along 2 : 
and along the spheres are k\\/ai and k^j 02 -, respectively. 

Then k\\, k 4 _^ 1, while a glance on Eqs. (3.3), (3.4), (3.5) and (3.9) implies 


L 





(5.1) 


where L is any linear combination with coefficients of order unity. Implementing these in the 
equations we get 


~ P ~ 




- 2 

ai 


(5.2) 
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and 


(5.3) 


D ~ C ~ —B 




’ 


£ = |c. 


Sc^-T.B 


3^2 \ k 


£ H—1 ^ 
2/02 


o? 


-1 


(5.4) 


Thus in the geometrical optics approximation the order of the dimensionless quantities relates 
as O (ajB) = O {afP) ^ O {C) = O {D) = O {aiE) = O {aiF) = O {ar^S) =0{l). 


5.1 High frequency evolutions of the uncoupled gravitational perturbations: 
gravitational waves 

The relevant coefficients are approximated as 


/kf, \ k?, 

9£0 - qm^{l-C)\ 4 + 7J-2 =-2 +-T’ 


~ 3(F+SC) - ^In 


qui ^ 3(F+SC) - 

at 



The damped wave equations (4.55)-(4.56) simplify to 


^fc||fcx+ + 


[ay al) 


= 0 , 


(5.5) 

(5.6) 

(5.7) 


(5.8) 



+ + 



(5.9) 


Both of these equations are of the form X + 2(T1X + VP'X = 0, where fl represents the 
undamped angular frequency. For C < 1 the oscillator is underdamped and the real angular 
frequency is given by The propagation speed of the wave therefore is Cw = 

— C'^/kphys (with kphys given by A:||/oi or k_i_/a 2 , when the propagation is along z or 
along the spheres). These considerations imply that we are assuming a negligible change in 
the scale factors over one period. 

In order to continue the analysis we define a small parameter e ~ {aikphys)~^ (char¬ 
acterising the geometrical optics approximation) and we consider perturbations along the 2 : 
direction and along the sphere separately. 


5.1.1 Waves propagating along the ^-direction 

We get C ~ 0 and T-|-SC'= ^ S, hence the propagation equations are 



20 3S 


ai 


X.T 

^kuk I 


kl 


+ 


a 


—T 

2^fc||fcx 


= 0 


xxT 


+ 


20 + 3S + — ) 
aij 


4/^ 


_i_ II 

'^kuk I 


= 0 


(5.10) 

(5.11) 
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The damping parameter turns out to be 


and the speed of propagation of both and is 

C|| = Vl - ^ 1 - Y = 1 - O {e^) ■ (5.13) 

~~T 

Thus to linear order in the geometrical optics approximation both ^^P^esent 

gravitational waves propagating with the speed of light. 


5.1.2 Waves propagating along the spheres 

We get C ~ —1 and F+T,C= hence 

” T f ” "X ' T' 

^ fc || fcx + (^20 + 0 ’ (^•^^) 

+ (^20 - Y + ^ 

Then there are different damping parameters for the two fields: 

and the speeds of propagation of different, but only to second order: 

C- 

c^-E^l-f = l-0{e^) , (5.18) 

c±w ^ l-^ = l-0(e") ■ (5.19) 


Again, to linear order in the geometrical optics approximation both £k\\k± ^®P" 

resent gravitational waves propagating with the speed of light. This is consistent with the 
generic theory, where the gravitational waves appear at the first order of the expansion of 
the Einstein equations. The second order terms, neglected in this picture can be interpreted 
as backreaction, leading to wavenumber-dependent dispersion. 


5.2 High frequency evolutions with matter sources: gravitational, shear and 
matter waves 


To first order in the geometrical optics approximation, the shorthand notations appearing 
in the evolution equations (4.71)-(4.73) simplify as follows. The coefficients of the algebraic 
terms of the gravitational perturbations are 


gso - qeo - Oho 




(5.20) 
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the coefficients of the damping terms become 


gsi ^ 2S + ^ , 

(^ + ^ + (2£: + 3S2) ^ ^ 

..1 = 0 + ^- 


, 2k}, u2 

s -4 + % 

n “2 


2^n 


ai I — Y “I —2 
\ «I «2 , 


^ (£ + ¥)f+ (2£ + 3S2)^ 

‘iu. = ^ + ^ - bb ,..x - "-s'" 


/Cn ^2 

II I ^_L 

“y “T —T 


ei^ + 4 

if al 


2^? fcf 

«i I T + ^ 


the coefficients in the algebraic source terms read 
seo = 


^ - I - (0 + 3S) - + A In [02 in + p)] 


seo = 


0 7S d 


1-c^ 
2 02 


^Wo = --—+—In 


dt 


V «I ^2 , 


while the coefficient of the time derivative source term in Eq. (4.72) is 

2^1? k} 


1 


sei = 


„2 „2 


-1 


d 

dt 


ai 


In ■ 


+ 

^1 0-2 


p + p 


-h 


(Zp 




2/Cm ^2 

s ^ + 4 

if a| 


Hence the evolutions (4.71)-(4.73), to leading order simplify as 


%\kl- + 


(1-c^) -V 

a2{pi+pf^^^^^ ’ 


ET 


k\\k 


-\-r/c-E'^ _L _L -n 


—T —T / /c| /j2 

The fourth equation of the closed system becomes 


= Q 


dknk, = a 2 id+P) 


all 


a? ^ o2 i ^11^^ 


2 ik 


(5.21) 

(5.22) 

(5.23) 

(5.24) 

(5.25) 

(5.26) 


(5.27) 

(5.28) 

(5.29) 

(5.30) 

(5.31) 
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Inserting Eq. (5.31) into Eq. (5.28) we obtain 


■‘knk. 


+ <lT,l^knk, + 



«i (1 - cj) 

2ik\\ 



(5.32) 


We comment on the system (5.29)-(5.32) as follows. The most striking feature is that the 
gravitational sector and fully decouples from the matter density gradient 

_ '-p 

While £^k^^ und Ti-k^^k^ obey damped oscillator equations (similarly to their counterparts 

p 

£k\\k± und discussed in the previous subsection), undergoes a forced oscillation. 

Let us first discuss the damped oscillations in the manner of the previous subsection. 


5.2.1 Gravitational waves 






+ 


1/2 


while the 


The undamped angular frequencies of £k^^k± und 'H^^kj_ £l = 

damping factors are (g = ^ and . For the waves propagating in the 2 and spherical 

directions, respectively, we get the following damping factors of O (e): 




2£ + 3E2 d ^ 


and 


C±£: — 
^±n ^ 


02 

2A:_l 

a 2 

2A:_l 


^ 2T + 3S2 d ni 

0 H-- In ^ , , 

2S dt alJ 

^ 2T + 3S2 d 


also the corresponding propagation speeds: 


1 - 




= 1-0(8^) , 


C||f ^ 

^ 1 - ^ = 1 - O (e") 


C^- 


and 


c±£ ^ 1 - ^ = 1 - O (e^) , 
% = l-0(e^) . 


^±H - 1 


(5.33) 

(5.34) 

(5.35) 

(5.36) 

(5.37) 

(5.38) 

(5.39) 


Thus, to leading order in the geometrical optics approximation the gravito-magnetic variables 
£k\^kj^ und represent pure gravitational waves. The non-identical corrections at higher 

order represent backreaction. 
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5.2.2 Shear waves and matter density gradient waves 

To leading order in the geometrical optics approximation Eq. (5.32) represents a wave for 

the shear propagating with the speed of sound Cg. At higher order both a damping 

mechanism and a force acts on this wave. The force is generated by the gravitational wave 
—T 

degree of freedom 

Again to leading order the shear is but a time derivative of the matter density 

gradient order to understand this claim it is necessary to remember that in the geo¬ 

metrical optics limit we consider wavelengths much shorter than the curvature radius, hence 
for the purpose of the wave propagation we can approximate the background scale factors 
and fluid characteristics as constants. Hence the latter also represents a wave propagating 
with the speed of sound Cg, but dephased with an angle 7r/2. This is in agreement with our 
assumption of an adiabatic speed of sound given by 

To higher order it mimics the damped and forced oscillation of . 


5.3 The degrees of freedom in the gravitational waves 


In the geometrical optics approximation we have obtained four equations representing grav¬ 
itational waves propagating with the speed of light, for the even and odd modes of the 2D 
electric and magnetic projections of the Weyl tensor, Sk^^k^i ■ Never¬ 

theless it is common knowledge that in general relativity gravitational waves carry only two 
degrees of freedom, represented by the -|- and x polarisations. In this subsection we address 
this apparent mismatch in the degree of freedom counting. 

We start by writing the geometrical optics limit of the uncoupled first order equations 
(4.48) and (4.49). By employing that at least one of the contitions S> 1 or A:_i_ ^ 1 holds 
together with O {a^B) = O {ajP) ^ O {C) = O {D) = O (aiE) = O (aiF) = O (a-^S) = 
0(1) and also the estimates (5.1), (5.2) and (5.3) we obtain 


ai [ay al) [al ^ al) 


(5.40) 



(5.41) 


Hence we have found that £k\ikj_ simply related, they represent the same 

degree of freedom. We have already shown that to leading order they obey undampened 
wave equations with the propagation speed of light. Again, in the geometrical optics limit 
the prefactors of the right hand sides can be considered constants, hence £k\ikj_ B-k^^kj^ 
are simply the time derivatives of each other, representing the same gravitational degree of 
freedom. Note that this analysis could have been done also for the general case in section 
4.3.1, where the same conclusions can be drawn from the system (4.48-4.49). However, since 
the conplings for the system in section 4.3.2 are more intricate, we have chosen to work in 
the geometrical optics limit throughout in this section. 

Next we revisit the geometrical optics limit of the coupled first order equations (4.63) 
and (4.64), which simplify as 



ai 

2 ik\ 




(5.42) 
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2ik\ 


-1 


(5.43) 



2ik\\ 


ai 





Thus, again, the prefactors on the right hand sides can be considered constants in the geomet¬ 
rical optics approximation, thus the Weyl variables and satisfying undampened 

wave equations with the propagation speed of light are the time derivatives of each other, 
representing the same gravitational degree of freedom. Note that in the geometrical optics 
limit this second gravitational degree of freedom also decoupled from matter and all four 
quantities T = obey 


3 > + 



T = 0 , 


(5.44) 


but they represent only two degrees of freedom. 

Like in the case of FLRW perturbations there are two matter degrees of freedom. One 
for the fluctuations of the density (via gradients of the energy conservation equation) and the 
other coming from the scalar part of the shear equation (representing velocity perturbations). 


6 Concluding Remarks 

A general treatment of vorticity-free, perfect fluid perturbations of Kantowski-Sachs models 
with a positive cosmological constant was considered within the framework of the 1+1+2 
covariant decomposition of spacetime. We showed that the system of pertnrbation equations 
can be organised into a hierarchy of three systems, namely (i) two conpled gravito-magnetic 
first order differential equations, (ii) four first order differential equations for the two com¬ 
plementary gravito-magnetic variables, a variable describing the shear of the world-lines and 
the gradient of the matter density perturbation, (hi) an extended set of algebraic relations 
involving all variables, which provides a way of determining their evolution. 

By assuming that the perturbation wavelengths is much smaller than the curvature 
radius of the Kantowski-Sachs background, we were able to use the geometrical optics ap¬ 
proximation to describe the evolution of high frequency pertnrbations. We found that system 
(i) gave rise to the leading order decoupled propagation equations for gravitational waves on 
this background, while to the next order, damping effects make the propagation along the 
spheres dephased. At leading order, system (ii) gives rise to two decoupled gravitational wave 
propagation equations for the complementary gravito-magnetic variables, supplemented by 
wavelike evolutions for both the shear and matter gradient pertnrbations, which both prop¬ 
agate with the same speed of sound < 1, out of phase by 7r/2. At the next order the 
gravito-magnetic oscillations are again damped, while the shear and matter waves obey forced 
oscillation wave equations. 

We note that the perfect fluid is marginally stable under the vorticity-free, anisotropic 
pressure-avoiding pertnrbations at high frequency. The degrees of freedom propagating as 
gravitational waves in the geometrical optics approximation are exactly the even and odd ten- 
sorial perturbations of both the electric and magnetic parts of the Weyl tensor, in agreement 
with its generic interpretation. While we have found four such quantities obeying undamp¬ 
ened wave equations with the propagation speed of light, the even electric and odd magnetic 
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Weyl projections represent the same gravitational degree of freedom, while 

the odd electric and even magnetic Weyl projections the other one. 

Beyond the geometrical optics approximation we have found indications for the existence 
of direction dependent dispersion relations. Remarkably, the second gravitational degree of 
freedom does not decouple from the matter density perturbation, unlike in 

Friedman universes. 
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A The relation between the 2D and 4D curvature tensors 


We give here the proof for the equivalency of Eqs. (2.26) and (2.27). The second covariant 
derivative of any 2D dual vector field 14 projected to the 2D subspace with gives 

= 6a6bVc + [{6^) (<5a«c) - (<5f>nz) {SaUc)] 

- (SaUb) 14 + {6anb) % . (A.l) 


Then from the definition of the Riemann tensor (VjVj — VjVj) 14 = RijkiV\ we find 

= 25^Jb]yc - 2 (<5[anfe]) 14 + 2 (<5[,nb]) Rg 

- [{^aUl) {6bUc) - (SaUc) {6bUl) 

- (Sani) {SbUc) + {SaUc) (^bU/)] R* . (A.2) 

With 6 [aUb] = ^£ab and 6 [anb] = the above identity reduces to 

5^ah]yc - ^eabV-c + iSaby-c = ^ [iV„W,'A/+ (5„Urf) (dftU,) 

- {SaUc) (SbUd) - {SaUd) {SbUc) + {SaUc) {SbUd)] ■ (A.3) 


The square bracket on the right hand side is the 2D curvature tensor TZabcd, as can be seen 
by comparing the left hand side with the dehnition (2.26). 


B Infinitesimal frame transformations on the Kantowski-Sachs background 
filled with perfect fiuid 

An infinitesimal frame transformation from the dyad (u“, n“) to the dyad (7t“, n“) can be 
defined as (see for higher dimensional spacetime [64]): 

Ua = Ua + Va + VUa , with = 0, (B.l) 

Ua = na + la + mUa , with u'^la = rfla = 0, (B.2) 

where Va, la, are all first order. We will neglect the second order contributions. The 

new dyad also obeys 

U^-Ua = -1 , n“na = 1 , Ifua = 0 , (B.3) 
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which implies 

V = m . (B.4) 

There are five gauge degrees of freedom to fix the frame on the perturbed spacetime. The 
transformations with Va = la = 0 represent 2D infinitesimal Lorentz boosts, while the pa¬ 
rameters Va and la are related to infinitesimal translations. 

The fundamental algebraic tensors Nab and £ab change accordingly: 


^ab ^ab T ‘^'^[a^b) ‘^^{alb) j 

^ab ^ab T ‘^'^[a^b]cl ‘^^[a^b\c^ 


(B.5) 

(B.6) 


The new 2-metric obeys Nabn^ = NabU^^ = 0. 

The kinematic quantities defined for the new dyad vectors arise from the decomposition 
of the covariant derivatives oiua and ria similarly to that given in Section 2. This implies the 
following transformations rules on Kantowski-Sachs background for the kinematic quantities: 


0 — Q V -\- 5aV°' , 

^ = -^ \ V 


1 


D — D -|- —£^^6aVb , 

'^ = + -I' — 5 


e 


20 


— *^a o ( ^ o ) ’ 


i^a = ^a- 
'^a = Sq — (Va + Sab' — 3S/a) , 


Oa — ®a T ^ T 


0 


(B.7) 

(B.8) 

(B.9) 

(B.IO) 

(B.ll) 

(B.12) 

(B.13) 

(B.14) 

(B.15) 

(B.16) 


— Ola T ia ; 

Sab = Sab + <5{a'Cb} > (B.18) 

Cab = Cab + S^ak} , (B.19) 

The gravito-electro-magnetic quantities £,71, £ab and Tiab are invariant under the in¬ 
finitesimal frame change, while the transformation laws of £a and Tia are 


7ia = 7ia- ‘^eabv'" 


(B.20) 

(B.21) 
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The matter variables /r, p, 11, Ha and Yiab are invariant under the infinitesimal frame 
change, while Q and Qa describing the energy current transform as 

Q = Q-{p + p)h', (B.22) 

Qa = Qa- {^J^■+P)Va . (B.23) 

The gauge-invariant variables defined by Eq. (4.1) transform as 

Ga = Ga + Gva , (B.24) 

where Ga = {pa, Xa,Va,Wa,Pa] and G = {/r, f, S, 0,p}, respectively. 


C Commutation relations 


The commutation relations of covariant derivatives of the scalar field T on Kantowski-Sachs 
background, to first order are 


4^ - T = ® ^ T , 

{5,^y = -A'k (s - , 

= -2eab^^i’ , 

Similar relations hold for the first order 2-vector Ta: 

® ^ Ts , 


(C.l) 

(C.2) 

(C.3) 

(C.4) 

(C.5) 


5ai>b - iV/iV/ (<5e^d)' = , (C.6) 

- IV/iV/= 0 , (C.7) 

5[A]^c = 2771Ve[a^6] , (C.8) 

and for the first order symmetric, trace-free 2-tensor Tq;,: 

^{ab} - ^{ab} = Ta , (C.9) 

•^a'kbc - N/N,^NJ (5,Te/)- = -^ (s - 6a^bc , (C.IO) 

Sa^bc - = 0 , (C.ll) 

25[a<Jb]^cd = 77 (A^c[a^b]d + 7Vd[a4'fe]c) ■ (C.12) 

where 77 is given by Eq. (3.9). 
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D Properties of vector and tensor spherical harmonics 


In this Appendix we enlist a set of identities for the even and odd Qa vector spherical 
harmonics, including the orthogonality relations 

= 0 , (D.l) 


the algebraic relations 

= -eM" , 

Qa = , 

(D.2) 

and the differential relations 

= 0 , 

—k± ^k± 

Qa =Qa = 0 , 

(D.3) 



(D.4) 

“2 

«2 


^aQk^ = 

02 


(D.5) 


= 0 , = —Q’"^ • (D-6) 

0.2 

The even and odd tensor spherical harmonics obey in turn the orthogonality relations 


the algebraic relations 




Q"at = - Q^at = > 


and the differential relations 


= Qlt = 0 , Qa:= Qa, = 0 , 


—kj_ 


= 


4 — 


-Qlt, ^"Qat = 


4 — 


= 


2-k 


-L/ofc 


202 


bT^^i 


Ql^ , 5°Qat = - 


2- k, —k 


202 


-Qa , 




f,7=rfcj 2- k"] 


202 




= 


2-k\—k 




202 


—fc, 2 - A;^ 


202 


-Q\ 


k± 


For reviews of various types of harmonics used in relativity see, e.g., [73, 74]. 


(D.7) 

(D.8) 

(D.9) 

(D.IO) 

(D.ll) 

(D.12) 

(D.13) 


E Harmonic expansion of the vorticity-free perturbation equations 

We give here the harmonic decomposition of Eqs. (4.3)-(4.34), employing Eq. (4.45) and 
= 0. The perturbation equations decouple into two sets, describing the even and odd 
parity sectors. 
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E.l Odd parity sector 


The evolution equations for gauge-invariant 2 -vector perturbation variables with odd parity 


are 


J.V 


iki 


ai 




/V 


0 


7T^ 
®fc||A: I 







kiik I 
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(E. 3 ) 


—V iku 

y — _ 

n^k^ 2 oi ^ 11 ^" 


+ T S- 




1 if 


38 
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(E. 4 ) 


The evolution equations for the odd parity tensor perturbations are 


■ T 



- y — 'T 
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(E. 5 ) 
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(E. 6 ) 

(E. 7 ) 

(E. 8 ) 

(E. 9 ) 

(E.IO) 

(E.ll) 

(E. 12 ) 

(E. 13 ) 

(E. 14 ) 

(E. 15 ) 

(E. 16 ) 
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E.2 Even parity sector 

The evolution equations governing the gauge-invariant scalar perturbation variables with 
even parity are 


9 k»kj 


24 " 


njS 


s - 




2Q\ 


0,5 _ _ ofrpS 


" ^ z, Y <J / " ^ U2 

The evolution of even parity 2-vector perturbations are given by 
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The evolution equations for the even parity 2-tensor perturbations 

yT -(y_^^yT ^ — Ay 
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The constraint equations for the even parity sector are 
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(E.30) 

(E.31) 

(E.32) 

(E.33) 

(E.34) 

(E.35) 
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We have checked the consistency of the evolution equations and constraints (Eqs. (E.35), 
(E.9) and (E.43), containing no evolutions for ■^\\k^ PX\\k^) showing that the 

time derivatives of the latter are identities. In fact not all algebraic relations are independent, 
their relations are summarised as follows: 1) Eq. (E.35) follows from Eqs. (E.9) and (E.43); 
2) Eq. (E.IO) follows from Eqs. (E.16), (E.33), (E.34) and (E.44); 3) Eq. (E.45) follows 
from Eqs. (E.12), (E.37), (E.41) and (E.46); 4) Eq. (E.40) follows from Eqs. (E.13), (E.15), 
(E.16) and (E.44); 5) Eq. (E.15) is consequence of Eqs. (E.ll), (E.14), (E.16), (E.33), (E.34) 
and (E.44); and finally 6) Eq. (E.32) is consequence of Eqs. (E.12), (E.37), (E.38), (E.39), 

‘■fciifcj.j 
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Imposing an equation of state p = p{p) gives ~ ^sPX^^k^ where = dp/dp is the 

adiabatic speed of sound squared. With the freedom in the choice of frame (Appendix B) 
the condition Oa = 0 can be set. In this case there are 19 constraints for 25 variables. By 
the frame choice and Eq. (E.33) we have = 0. The evolution of all other harmonic 
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where B and C are defined in Eqs. (4.50) and (4.51), respectively, and 
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